APPARATUS AND METHODS FOR 
DIAGNOSING OSTEOPOROSIS AND OTHER DISEASES 
WITH MR IMAGING 

FIELD OF THE INVENTION 

This invention generally relates to magnetic resonance imaging 
(MRI) utilizing nuclear magnetic resonance (NMR) phenomena. It more 
particularly relates to a specialized MRI system used for diagnosis 
of osteoporosis and other diseases. 

BACKGROUND OF THE INVENTION 

Osteoporosis is a widespread disease characterized by bone loss 
usually associated with aging. It is a complex and chronic disease 
that may not be recognized until resulting fractures late in life. 
In many countries around 50% of women above age of 70 is osteoporosis 
and to overcome it has become an important social requirement. 

For diagnosis and follow-up of osteoporosis, measurement of bone 
mineral density (BMD) is indispensable. The conventional methods 
for the BMD measurement are based on X-ray attenuation measurements 
through a bone. Among them, DEXA (Dual-Energy X-ray Absorptiometry) 
is now the standard method, but cannot be used repeatedly nor for 
young women because of the ionizing radiation. Thus the X-ray BMD 
measurement should be replaced by some noninvasive and high precision 
method such as bone marrow volume measurement using MRI as described 
below, because MRI has no ionizing radiation and spatial localization 
capability is guaranteed. 

The bones of a human skeletal system have a dense outer shell made 



of cortical bone (alternatively termed "compact" or "dense" bone). 
Inside the cortical bone, many regions have a mesh or network of 
trabecular bone (alternatively termed "spongy" or "cancellous" bone) , 
made up of roughly the same material as the cortical bone. The regions 
enclosed by the cortical bone and/or trabecular network are filled 
with "bone marrow", of which major components are water and lipid. 
Since the protons (hydrogen nuclei) of the water and lipid can be 
measured using MRI , the bone marrow volume can be quantified using 
MR images . Thus quantitative bone MRI can be a complementary method 
for the measurement of BMD because bone is almost chemically 
homogeneous. Studies according to this concept are found in the 
following references: Jara el al. , "High-Resolution Variable Flip 
Angle 3D MR Imaging of Trabecular Micros true ture in Vivo" , Magnetic 
Resonance in Medicine 29: 528-539, 1993; Wehrli et al., "Cancellous 
Bone Volume and Structure in the Forearm: Noninvasive Assessment 
with MR Microimaging and Image Processing", Radiology 206: 347- 
357, 1998; L. Hilaire et al. , "High-Speed Spectroscopic Imaging for 
Cancellous Bone Marrow R 2 * mapping and Lipid Quantification", 
Magnetic Resonance Imaging, 18: 777-786, 2000. 

MRI is by now a commercially available and widely accepted 
non-invasive method for measuring information about internal 
structures of living systems and heterogeneous materials. However, 
most of medical MRI are developed for a whole human body and require 
a large cost and space. For the measurement of bone marrow volume, 
a compact MRI system for a part of a human body is appropriate because 
bone density for a whole human body can be well estimated from a 
measurement of a part of bone. 
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SUMMARY OF THE INVENTION 

The most serious fracture for osteoporosis patients is the hip 
> fracture because the patients cannot move by themselves for more 
than several months. To acquire MR images of a hip bone (femur) for 
diagnosis of osteoporosis, a whole body MRI system is indispensable, 
requiring a large cost for the patients. However, because 
correlations among densities of bones at various parts of a human 
jj 10 body are very high , the bone density at the femur can be well estimated 
from those of other bones. 

In this invention, calcaneus is used for the bone density measurement 
using MRI because of the following reasons: first, calcaneus is made 
of homogeneous trabecular bone and reflects the whole -body bone 
''4 15 status very well; second, calcaneus is a peripheral part of a human 
3 body and easy to measure using a compact MRI system; third, the magnet , 

*** usually the heaviest unit of an MRI system, can be placed in the 

lower level and the MRI system becomes mechanically stable. For the 
measurement of a heel, an RF probe box having an oval aperture with 
20 a support pad to fix the heel in the RF coil is used. In the RF coil, 
external reference materials are placed and imaged together with 
the heel for NMR signal intensity calibration. 

Quantification of materials using a compact MRI system gives rise 
to several new problems, because the magnet, gradient coil, and RF 
25 coil are small compared to those of a whole body MRI and the 
inhomogeneity of the static magnetic field, magnetic field gradients , 
and RF magnetic field become large compared to that of a whole body 
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MRI . To measure proton density of bone marrow of calcaneus in such 
inhomogeneous magnetic fields, spin-echo imaging sequences with 
different echo times are used to image a heel together with external 
reference materials. By combining measurements with different echo 
5 times and reference materials , T 2 decay effects and spatially varying 
image intensities are corrected. As a result, the proton density 
of bone marrow is obtained and used for evaluation and follow-up 
of osteoporosis and other diseases. 

Jj 10 BRIEF DESCRIPTIONS OF THE DRAWINGS 

|1 FIG.l is a schematic overview of a compact MRI system for imaging 
of a heel; 

ss'jst 

m FIG .2 is a schematic representation of an RF probe box with an oval 

^ aperture used for imaging of a heel; 

^ 15 FIG. 3 schematically depicts a solenoid coil and matching circuit 

£3 for the RF probe; 
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FIG. 4 is a schematic cross -sectional depiction of a support pad for 
a heel and external reference materials, inserted in the open-bore 
of the RF probe shown in FIG. 2; 

FIG. 5 is a schematic cross -sectional depiction of another type of 
the support pad for a heel and external reference material, inserted 
in the open -bore of the RF probe shown in FIG. 2; 

FIG. 6 is a standard 2D spin-echo pulse sequence used to quantify 
the proton density of bone marrow in calcaneus; 

FIG. 7 is a flow diagram of one method for calculating proton density 
using the MRI system shown in FIG.l; 

FIG. 8 schematically shows a vertical cross -sectional image of a heel 
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with external reference phantoms; 

FIG. 9 schematically shows a vertical cross -sectional image of a water 
phantom for intensity correction of images of heels with external 
reference phantoms; 
5 FIG. 10 is a semilogarithmic plot of spin-echo image intensities of 
a plant oil phantom plotted against the echo time; 
FIG. 11 is a large flip-angle excitation 2D spin-echo pulse sequence 
to shorten the repetition time for the quantification of the proton 
density of bone marrow in calcaneus. 

O io 

^3 DETAILED DESCRIPTION OF THE INVENTION 
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FIG.l schematically illustrates an MRI system with a permanent 
magnet, gradient coil, RF probe, and MRI console for imaging of a 
15 heel, as an exemplary embodiment of this invention. 

In FIG.l, the heel 1 is inserted into the RF coil 2 placed between 
fc* the permanent magnet 3 and gradient coil 4. The permanent magnet 

produces a homogeneous magnetic field of about 0.2 IT over the 12 
cm diameter spherical volume in the central region of the magnet 
20 gap, of which width is 16 cm. The volume of the homogeneous magnetic 
field is sufficient to image a heel and external reference phantoms 
at the same time. The RF coil is actually stored in an RF shielded 
box as shown in FIG. 2 to eliminate the electromagnetic coupling 
between the RF coil 2 and gradient coil 4. 
25 MRI measurements are controlled by the MRI console, which consists 
of a display monitor 5, computer system 6, RF transceiver 7, gradient 
driver 8 , and RF transmitter 9 . The computer system includes a digital 
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signal processor (DSP) board for the MRI pulse programmer, a direct 
digital synthesizer (DDS) board for generation of a radio frequency 
source for the NMR transceiver, and an analog -to- digital converter 
(ADC) board for digitization of the NMR signal. 

The RF transceiver 7 is used to generate RF pulses for excitation 
of the proton spins of the heel and convert the NMR signal modulated 
at the Larmor frequency at around 8 . 9 MHz to a lower frequency ( several 
tens of kHz) having spatial information of the proton spins. The 
gradient driver 8 is used to supply electric currents to the gradient 
coil 4 according to waveforms generated by the pulse programmer. 
The transmitter 9 is used to supply RF currents to the RF coil 2. 

FIG. 2 schematically illustrate the RF probe used in the system. The 
RF probe box 10 is made of brass plates with an oval aperture 11 
in the upper side. The dimension of the box is 22 cm (width) x 16 
cm (height) x 13 cm (depth) and that of the aperture is 16 cm (long 
axis ) x 8 cm ( short axis ) . The aperture can accommodate heels of 
most female adult patients as depicted as 12. The probe has an 
open-bore RF coil support 13 made of acrylic plates, which goes from 
the upper aperture down to another aperture in the lower side of 
the probe box. A solenoid coil with an oval cross -section is wound 
over the coil support. In the open-bore coil support, a support pad 
for a heel and external reference phantoms is inserted as will be 
depicted in FIG. 4 and FIG. 5. An RF connector 14 is attached on the 
upper side of the probe box to supply RF currents for excitation 
of proton spins and to receive NMR signals from proton spins. 

An LC tank circuit is schematically depicted in FIG. 3 . In this figure , 
although the RF coil 15 is depicted as a solenoid with a circular 



cross -section, the solenoid coil with the oval cross-section which 
is wound on the RF coil support 13 can be also used for the tank 
circuit . The capacitor 16 is used mainly for tuning of the tank circuit 
to the resonance frequency (8.9 MHz) and the capacitor 17 is used 
» mainly for impedance matching to 50 ohms. 

FIG. 4 is a schematic cross -sectional depiction of a support pad 18 
for a heel 19 and external reference materials denoted by 20 and 
21, inserted in the open-bore of the RF probe shown in FIG. 2. For 
quantification of proton densities of bone marrow in calcaneus, to 
P 10 fix the heel and to acquire MRI images together with known reference 
^ materials such as water is essential, because tuning of the tank 

Sft circuit varies person to person dependent on the size and position 

of the foot. In this figure, polystyrene foam is used for the support 

if\ 

3 pad 18 and CuS0 4 water solution in cylindrical polyethylene bottles 

C3 

Nl 15 is used as the external reference material. 

FIG. 5 shows an alternative embodiment for the support pad 22. A 
heel 23 is directly placed on the external reference material 24 
which is flexible to fit any shape of a foot to improve the precision 
of RF field inhomogeneity correction. 
20 FIG. 6 is a standard 2D spin-echo pulse sequence used to quantify 
the proton density of bone marrow in calcaneus . Image intensity I(x,y) 
acquired with the spin- echo sequence can be expressed except the 
J modulation effect of lipid as follows: 

25 I(x,y) = kf(x,y)p(x,y){l-p(x,y)zxp(-TR VT^y^joxp^TE /T 2 (x 7 y)) , (1) 

where k is a constant, f(x,y) represents spatial variation of image 
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intensity for a uniform sample, which is determined by inhomogeneity 
of static magnetic field, magnetic field gradients , RF magnetic field, 
and coil sensitivity (in principle the same as distribution of RF 
magnetic field) , p(x,y) represents a factor reflecting longitudinal 
5 magnetization, which is unity for perfect 90 degree excitation over 
the cross- sectional plane, 77?* is a time close to the repetition 
time TR, T x {x,y) and T 2 (x,y) are T x and T 2 distributions in the plane. 
FIG. 7 shows a flow diagram of one method to quantify the proton 
density of bone marrow in calcaneus . To quantify the proton density 
M 10 p(x,y) / spin-echo sequences with TR/TE= 1200ms /12ms and 



TR/TE=1200ms/108ms are used for the heel measurements as denoted 
by 25. Since the T x of the protons of the bone marrow in calaneus 

Kip 

<j£ is about 250 ms at 8 . 9 MHz , MR image intensities at calcaneus acquired 

with the above sequences can be expressed as 
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I l (x J y)=kf (x, y ) p(x, y ) exp(-12 /T 2 (x,y)) ( 2 ) 

I 2 (x,y) = ty(x,y)p(x 9 y)exp(-108/T 2 (x,y)) . (3) 



20 By using above two equations, T 2 (x 9 y) and kf (x,y)p(x,y) are calculated 
as denoted by 26 and 27. For the measurement of kf(x,y) , a CuS0 4 water 
solution phantom of which T 2 is about 43ms is placed in the RF coil 
instead of a heel and imaged with the same spin-echo sequences with 
TR/TE=500ms/12ms and TR/TE=500ms/48ms as denoted by 28. By using 

25 these images, kf(x,y) is calculated when the proton density of water 
is defined as unity as denoted by 29 and 30. The proton density p(x,y) 
of the bone marrow when the proton density of water is defined as 
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unity is thus calculated by dividing kf (x,y)p(x 9 y) bykf(x,y) as denoted 
by 31. The trabecular bone volume fraction is then computed as 
1- p(x,y) , as denoted by 32. 
FIG. 8 and FIG. 9 schematically show cross -sectional images of a heel 
5 and a water phantom with external reference phantoms. For actual 
calculations for proton densities and T 2 relaxation times, mean 
values in the squares depicted in FIG. 8 and FIG. 9 are used because 
random noise is superimposed on the MR images and the precision is 
limited . 

S3 10 FIG. 10 is a semilogarithmic plot of spin-echo image intensities of 

%y 

feS a plant oil phantom measured against the echo time, showing the J 

iP modulation effect on the spin-echo intensities . The similar property 

# is observed for the bone marrow in calcaneus because the bone marrow 

# has a similar chemical composition as the plant oil. The J modulation 
%j 15 effect, which is seen as the deviation of the spin-echo signal shown 

w 

□ from the straight line, is shown as an oscillation with about 100ms 

y* period. Thus the T 2 decay of the lipid protons can be corrected 

reasonably using the images acquired with TE=12ms and TE=108ms. 
FIG. 11 is a large flip-angle excitation 2D spin-echo pulse sequence 
20 to shorten the repetition time for the quantification of the proton 
density of bone marrow in calcaneus. By setting the excitation flip 
angle to an angle between 90 and 180 degree, the recover of the 
longitudinal magnetization becomes faster as depicted in FIG. 11. 
The arrows denoted by 33 to 41 show a time variation of longitudinal 
25 nuclear magnetization. 

While only a few specific exemplary embodiments of this invention 
have been described in detail, those skilled in the art will readily 
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appreciate that many variations and modifications may be made in 
these exemplary embodiments while yet retaining many of the novel 
features and advantages of this invention. Accordingly, all such 
modifications and variations are intended to be included within the 
scope of the appended claims . 
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